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Abstract. An anomalous nuclear modification was reported by JLab measurements on the 
beryllium-9 structure function F2. It is unexpected in the sense that a nuclear modification slope 
is too large to be expected from its average nuclear density. We investigated whether it is explained 
by a nuclear clustering configuration in ^Be with two a nuclei and surrounding neutron clouds. 
Such clustering aspects are studied by using antisymmetrized molecular dynamics (AMD) and also 
by a simple shell model for comparison. We consider that nuclear structure functions F2 consist of 
a mean conventional part and a remaining one depending on the maximum local density. The first 
mean part does not show a significant cluster effect on F2. However, we propose that the remaining 
one could explain the anonymous JLab slope, and it is associated with high densities created by the 
cluster formation in ^Be. The JLab measurement is possibly the first signature of clustering effects 
in high-energy nuclear reactions. A responsible physics could be an internal nucleon modification, 
which is caused by the high densities due to the cluster configuration. 
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INTRODUCTION 

Nuclear modifications of structure functions F2 were found by the European muon 
collaboration, and it is often called the EMC effect. The nuclear modifications are now 
experimentally measured from relatively small x to large x. Theoretical mechanisms 
are different depending on the x region [1]. At small x, the modifications are caused 
by nuclear shadowing due to multiple scattering of a qq pair coming from the virtual 
photon. At medium x, nuclear effects are due to binding and a possible internal nucleon 
modification. Nucleon Fermi motion and short-range nucleon-nucleon correlation cause 
large-x nuclear effects. Using nuclear F2 data together with the other ones, optimum 
parton distribution functions have been determined in nuclei [2, 3]. 

Nuclear modifications at x > 0.2 are generally described by a convolution model with 
the nucleon structure function F2 convoluted with the nucleon momentum distribution 
in a nucleus. The nucleon four-momentum distribution, which is called the spectral 
function, has been calculated in a conventional shell model possibly with the short- 
range correlations and internal nucleon modifications. On the other hand, it is known 
in low-energy nuclear physics that some nuclei have clustering configurations which 
cannot be described by simple shell models. It is an an interesting topic to investigate 
such clustering aspects in deep-inelastic structure functions [4]. 
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FIGURE 1. Nuclear modification 

slopes \d{F^/F^)/dx\ by JLab. 



An anomalous nuclear modification was found for 
the beryllium-9 nucleus at the Thomas Jefferson Na- 
tional Accelerator Facility (JLab) by measuring ra- 
tios , where A and D denote a nucleus and the 
deuteron. Usually, nuclear modifications are smooth 
functions of the average nuclear density. JLab mea- 
surements showed the nuclear modification slope 
d{F2 1^2)1 dx by approximating their data by 
straight lines at 0.35 <x< 0.7 [5]. They are shown in 
Fig. 1 as a function of the scaled nuclear density. It is 
obvious from the figure that the nuclear modification 
slope of ^Be is too large to be expected from the average density. 

In Ref. [4], we proposed that the anomalous JLab result could be interpreted by a 
clustering phenomenon in the ^Be nucleus. Here, the clustering means that the ^Be 
consists of two a (^He) nuclei with extra-neutron clouds. Such a nuclear clustering 
phenomenon had never been investigated in the structure functions, although small quark 
clusters such as a six-quark bag model was considered in the early stage of EMC-effect 
studies. We explain that the JLab measurement could be interpreted if the modification 
slope is plotted as a function of the maximum local density at a cluster by using a 
theoretical method of antisymmetrized molecular dynamics (AMD). In this article, we 
explain our theoretical approach for explaining the anomalous JLab result. 




NUCLEAR CLUSTERING EFFECTS IN F2 

Nuclear structure functions at x > 0.2 are de- 
scribed by a convolution description shown in Fig. 
2. The nuclear structure function F^ is theoreti- 
cally evaluated in two steps: first by calculating 
nucleon {N) momentum distribution in a nucleus 
(A), second by calculating the quark {q) momen- 
tum distribution in the nucleon. Then, the virtual 

photon from the charged lepton (electron in the ^^r-^.^^ ^ ^ w ^ • ^ t-a 
. ^ . .11 1 1 4 FIGURE 2. Convolution description of i^f. 

JLab case) mteracts with the quark. Namely, F^ 

is given by the nucleonic one F^ convoluted with the nucleon momentum distribution 
in the nucleus f{y) [1,6]: 

F,\x,Q') = iyyf{y)F,^{x/y,Q^), M = \ j PNy 5 {y (1) 

where y is the momentum fraction y — Ma PN'^il {^N PA'4)' 

In order to illustrate clustering effects on F2 , we calculate the momentum density 
\^{pn)\^ in two methods, a simple shell model with the harmonic-oscillator potential 
M^^co^r^ /2 and the AMD. It is known that the wave functions of the harmonic-oscillator 
potential are expressed by Laguerre polynomials and spherical harmonics. The AMD 
(or FMD (fermionic molecular dynamics)) is a typical method in investigating cluster- 
ing aspects of nuclear structure. The advantage of the AMD is that it does not assume 



any specific structure, cluster- or shell-like configuration, on nuclei. The AMD is a vari- 
ational method, in which a nuclear wave function is given by the Slater determinant of 
single-particle wave packets: |0(ri,r2, • • -^ta)) = det[(pi(ri), (jC)2(r2), • • (Pa(^a)]/VA!. 
Here, (Pi{rj) is the single-particle wave function given by (Pi{rj) = (l>i{7j) Xi'^i with 
spin and isospin states Xi ^i^d T/. The function 0/(ry) is the space part expressed as 
0/(ry) = (2v / ;r)^/^ txp[—v{fj — %j V^)^] • Supplying simple A^A^ interactions, we deter- 
mine the variational parameters v and % by minimizing the system energy. Calculated 
AMD space densities are shown in Figs. 3 and 4 for ^He and ^Be, respectively. The 
^He AMD density is almost the same as the shell-model one. However, as shown in Fig. 
4, the ^Be density is totally different from a monotonic shell-model density because the 
^Be has the configuration of two a clusters with surrounding neutron clouds. Using both 
densities of the shell and AMD models, we calculate in Eq. (1). 




FIGURE 3. Density distribution of "^He [4] . FIGURE 4. Density distribution of ^Be by AMD [4] . 



Although the two-dimensional (x, y) den- 
sity for ^Be in Fig. 4 is very different from 
a shell-model one, inhomogeneity of the nu- 
clear density is washed out if the angular aver- 
age is taken [4]. Then, transforming the den- 
sities to the momentum space, we obtain the 
results in Fig. 5 for "^He and ^Be. It is inter- 
esting to find that the AMD wave function of 
^Be has larger high-momentum components 
than the shell-model one, whereas both densi- 
ties are the same in ^He. The high-momentum 
components are created by the cluster devel- 
opment in ^Be because nucleons are confined 
mainly in the two small-space clusters. 

Using the momentum densities, we cal- 
culate the nuclear structure functions at 
2^=5 GeV^ by using Eq.(l). The results are 
shown in Fig. 6 for Fr^^ jF^ together with ex- 
perimental data. Since the short-range corre- 
lations and internal nucleon modifications are 
not taken into account in our formalism, the 
curves do not fully agree with the data. How- 




FIGURE 5. Nucleon-momentum distributions 
in ^He and ^Be by shell and AMD models [4]. 




FIGURE 6. Nuclear modifications in ^Be 
by shell and AMD models [4]. 



ever, it is interesting to find some clustering effects in ^Be as shown in Fig. 6 because 
there are differences between the two curves of the shell and AMD models, whereas 
both are the same in ^He. If the current experimental errors and other effects such as 
the correlations are considered, it may not be easy to find a clear clustering signature in 
^1 1^1 fr^^ ^he conventional mean part. 

We found that the conventional binding and 




Fermi-motion contributions do not show a sig- ^ jLab(2oo9) ''c(amd) 

nificant clustering effect. Next, we consider the - ^ ^ ^; ' ^ 

possibility that the high densities due to the clus- ^ _ 

^ ^ i AMD 

ter formation in ^Be could be the origin for the | sheii 
anomalous JLab data. We plotted the JLab data by 
the maximum local densities calculated by both 

models in Fig. 7. The curve interpolates the shell- * 0.05 0.1 0.15 0.2 ^25 

model data points except for ^Be. If the ^Be slope Maximum locai density: maxfpw] a/W) 

data is shown by the maximum density of the shell FIGURE 7. Nuclear modification slopes 
model, it is again too large to be expected from shown by maximum local densities [4]. 
other nuclei. However, if it is shown by the max- 
imum density of the AMD with the cluster structure, it is on the curve. In this way, 
the "anomalous" JLab result could be explained by the cluster structure in the nucleus. 
There is also a small difference between the AMD and shell model in ^^C of Fig. 7, 
which is caused by the mixing of a cluster-like configuration in the AMD model of ^^C. 

Since it may be confusing for the reader that the ^Be slope is understood by the cluster 
structure although the effect is rather small in Fig. 6, we would like to explain our 
viewpoint. The nuclear structure functions consist of the mean conventional part and 
the remaining one depending on the maximum local density: 

F2 — (mean part) + (part created by large densities due to cluster formation) . (2) 

The first part is described by the usual convolution calculation with the spectral function 
given by the averaged nuclear density distribution. The remaining part is associated with 
the inhomogeneity of the nuclear density, before taking the average of nuclear wave 
function, given by the nuclear cluster structure. The physics behind the latter part could 
be a nuclear-medium modification of internal nucleon structure. Our studies suggest that 
the physics mechanism, associated with the high densities created by the clusters in ^Be, 
could be the origin for the anomalous slope d{Fr^^ jF^^jdx, Such studies of possible 
cluster structure in deep inelastic scattering will be continued at JLab [7], and a more 
elaborated theoretical model needs to be developed for comparison with future data. 
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